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We performed a systematic density functional study of the adsorption of copper, silver, and gold 
adatoms on graphene, especially accounting for van der Waals interactions by the vdW-DF and 
the PBE+D2 methods. In particular, we analyze the preferred adsorption site (among top, bridge, 
and hollow positions) together with the corresponding distortion of the graphene sheet and identify 
diffusion paths. Both vdW schemes show that the coinage metal atoms do bind to the graphene sheet 
and that in some cases the buckling of the graphene can be significant. The results for silver are at 
variance with those obtained with GGA, which gives no binding in this case. However, we observe 
some quantitative differences between the vdW-DF and the PBE+D2 methods. For instance the 
adsorption energies calculated with the PBE+D2 method are systematically higher than the ones 
obtained with vdW-DF. Moreover, the equilibrium distances computed with PBE+D2 are shorter 
than those calculated with the vdW-DF method. 

PACS numbers: 68.43.Fg, 71.15.Mb, 73.22.-f 



I. INTRODUCTION 

As of today, it is still a challenge for theory to correctly 
predict the binding of atoms and molecules on coinage 
metal surfaces!^ Recently reported studies, employing a 
variety of methods, showed the importance of non-local 
correlations for the correct description of the physisorp- 
tion of aromatic molecules on copper, silver, and gold 
surfaces Also the reverse setup of coinage metal atoms 
and clusters on carbon-based materials such as graphite 
is a long-standing and not completely resolved issue. Al- 
ready Darby et al. [H] speculated that dispersion forces, 
i.e. van der Waals forces, significantly contribute to the 
adsorption energies of gold atoms on graphite. 

In recent years graphene, the two-dimensional building 
block of g raphite, has gained much attention in its own 
right pHEO The gapless ultrarelativistic energy spectrum of 
graphene^ as well as effects such as Klein tunneling^ are 
of fundamental interest from a theoretical point of view. 
Moreover, graphene shows the Quantum Hall effect at 
room temperature,^ and its conductive properties cor- 
respond to ballistic transport. Also, graphene offers enor- 
mous possibilities for applications in elec tronic s, sensors, 
biodevices, catalysis, energy storage etc ! 11 * 14 ! that often 
involve adsorbed coinage metal clusters. For the utiliza- 
tion of these systems a deeper understanding of their sta- 
bility and electronic properties is needed. Adatoms and 
small clusters of coinage metals are known from experi- 
ment to be weakly bound and highly mobile on graphite 
and graphene JEHU] 

A wide range of theoretical studies have been per- 
formed on small coinage metal clusters on graphite and 
graphene ! 22 " 28 ! We are not aware of any experimental 
study of the adsorption sites or energies of coinage metal 
adatoms on graphene. However, there exists, to our 



knowledge, only one theoretical study that attempted 
to evaluate the importance of non-local correlation ef- 
fects on the binding of silver and gold nanoclusters on 
graphiteP^ The authors of Ref. [29] used a number of 
different semi-empirical methods, based on pair poten- 
tials, and concluded that the obtained results are in fact 
strongly dependent on the choice of the method, calling 
for more studies on these systems. 



To fill this gap, we undertook a systematic investiga- 
tion of the adsorption of copper, silver, and gold adatoms 
on graphene, based on density functional theory (DFT) 
and accounting for van der Waals (vdW) interactions by 
the vdW-DF method! ^HH as well as by the PBE+D2 
method^ In principle, more advanced schemes such as 
the Random Phase Approximation^ ^ (RP A) or the 
Quantum Monte-Carlo (QMC) methocP^ could also 
be used to take the van der Waals interaction into ac- 
count. Recently, the vdW-DF method has been com- 
pared to the RPA for the case of graphiteP^ It was found 
that, although the agreement was not perfect, the vdW- 
DF appeared as a quantitative method to take disper- 
sive interactions into account. The same comparison can 
be made between PBE+D2 and RPA for the same sys- 
tem. There, the binding energy for graphite was found to 
be 55 meV/atom with the PBE+D2,E2lto compare with 
the 48 meV/atom obtained with the RPA. Although the 
PBE+D2 and vdW-DF methods do not reach the accu- 
racy of advanced quantum chemistry methods, they are 
at the moment the only computational tools that allows 
to treat the dispersive interaction in such a large system 
as a sheet of graphene with a metallic adatom, together 
with periodic boundary conditions. 
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II. METHODS 

We calculated the total energies of Cu, Ag, and Au 
on either of the top (t), bridge (b), or hollow (h) po- 
sition on a 5x5 graphene sheet, i.e. 50 carbon atoms, 
using the calculated C-C bond length of 1.42 A. While 
the positions of the adatoms were fixed during all calcu- 
lations, the carbon atoms in the graphene sheet, except 
those in the rim of the supercell, were free to relax. The 
relaxation procedure was stopped, when the Hcllmann- 
Feynman forces on the carbon atoms, that were allowed 
to relaxed, were smaller than 5 • 10 -3 eV/A. 

These scalar-relativistic ab-initio DFT calculations 
were perfor med using the projector augme nted wave 
(PAWjSH m ethod as implemented in VASPpEH The 
exchange-correlation interaction was treated in the gen- 
eralized gradient approximation (GGA) in the param- 
eterization of Perdew, Burke, and Ernzerhof (PBE) 45 
and, for comparison, also in the local density approx- 
imation (LDA) in the parametrization of Perdew and 
ZungerpS A cut-off energy of 600 eV was used and a 
Gaussian smearing with a width of a = 0.05 eV for the 
occupation of the electronic levels. Spin-polarization was 
taken into account for all the calculations. The repeated 
graphene sheets were separated from each other by 20 A 
of vacuum. A Monkhorst-Pack T-centered 5 x 5 x 1 k- 
point mesh (13 k-points in the irreducible wedge of the 
Brillouin-Zone) was used for the structural relaxations of 
the carbon atoms. We also tested a finer k-point mesh, 
i.e. 16 x 16 x 1 k-point mesh (130 k-points), and found 
the changes in the geometry and the charge density dis- 
tribution of the systems to be negligible. 

We accounted for the non-local correlation energies by 
employing two different methods: the van der Waals den- 
sity functional (vdW-DF) methocPSHU! as implemented 
in the Jiilich Non-Local (JuNoLo) codep^ and the den- 
sity functional theory plus long-range dispersion correc- 
tion (DFT+D2) method^ in the implementation of Ref. 
[48) . From a practical computational point of view, these 
methods differ significantly. While vdW-DF is imple- 
mented as an ab-initio post-processing method in the 
JuNoLo code, the PBE+D2 method allows to include 
the vdW interactions in the self-consistency cycle, i.e. 
during structural relaxations. In Ref. (49) it has been 
shown, for instance, for a set of small molecular systems 
that the total energies and binding distances obtained 
with the vdW-DF as a post-processing method and in 
a self-consistent implementation agree very well. As an 
input to the vdW-DF calculations one uses the charge 
density of the relaxed structure obtained in DFT with 
PBE. In PBE+D2, as well as in PBE, we started from 
the same initial configuration, i.e. a flat graphene sheet 
with an adatom fixed above one of the adsorption sites. 
During the structural relaxation, all the carbon atoms, 
except for those in the rim, were free to relax. Note that 
in PBE+D2 the van der Waals forces were taken into 
account during the relaxation. 

Within the vdW-DF method, the missing vdW inter- 



actions are accounted for by replacing the correlation en- 
ergy from the GGA calculation by a sum of a strictly 
local E^ BA and a non-local correlation E™ 1 energy: 

^.vdW-DF _ ^PBE _ ^;PBE _j_ jtjPBE _|_ ^LDA _|_ ^nl ^ 

where Eq BE and £^ C BE are the total energy and the 
exchange-correlation energy, respectively, calculated with 
PBE. The non-local contribution to the correlation en- 
ergy is calculated as: 



E? = -J d 3 rd 3 r'n(f)(t)(r,f')n(r'), (2) 

with n(r) being the charge density. The kernel 4>(f,r') 
depends on the distance (r— r ') and the electron densities 
n in the vicinity r and r ' . 

In the PBE+D2 method, on the other hand, the 
vdW interactions are described by a pair-wise correc- 
tion (-Edisp)j optimized for some popular DFT function- 
als, and added to the self-consistent Kohn Sham energy 
(-Eks-dft) such that: 

£pBE+D2 = i?KS-DFT + Edi sp . (3) 

with 



"at "at ^ij 

= -fEEE V° -Vi° /(|r<, ° - rhL ^ 

»=i i=i l 1 1 

(4) 

where the summations are over all the atoms and all the 
translations of the unit cell. The prime indicates that 
i =/= j for L=0, ensuring that there is no double counting. 
The scaling factor sg is dependent on the Kohn-Sham 
functional. Its value is 0.75 for the PBE functional, which 
was used here. Cg J is the vdW coefficient for an atom pair 
(ij) and r l > L is the position vector of atom i after per- 
forming L translations of the unit cell along lattice vec- 
tors. Also, f(r) is a damping function which cancels the 
contribution from -Edisp f° r distances corresponding to 
standard (covalent or ionic) bonds. In practice, the sum- 
mations over all the atoms are replaced by summations 
up to a suitably chosen radius. All the pair interactions 
up to a radius of 12 A were included in the calculation of 
i?disp- Moreover, since gold was not included in the list of 
the elements provided in the original DFT+D2 article^, 
we have used a value of 40.62 Jnm 6 /mol for the Cq coef- 
ficient and of 1.772 A for the vdW radius of AuPl Note 
that for the DFT part (-Eks-dft) of the calculation, we 
used the same scalar-relativistic DFT method and PBE 
functional as described above. 

The adsorption energies -E a d s of a metal atom M on a 
graphene sheet G are calculated as: 
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A. Cu on graphene 
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where the Eo are the ground state energies of the adatom 
on graphene [M/G], graphene [G] and metal [M] alone, 
calculated with the marked method. The -E a ds are nega- 
tive when the adsorption is exothermic. 



III. RESULTS 

In this section, we present our results concerning the 
adsorption energies, preferred sites, and equilibrium dis- 
tances for Cu, Ag, and Au adatoms on graphene. In 
order to reach a qualitative conclusion, we used two dif- 
ferent schemes to take the dispersive interaction into ac- 
count. Namely, the vdW-DF and the PBE+D2 methods 
were employed, see the previous section for details. For 
the sake of comparison, we also report our results ob- 
tained with a semi-local (GGA) and a local (LDA) ap- 
proximation, although none of these methods are able to 
describe the non-local correlation effects needed to tackle 
the present problem correctly. 

The total energy curves, where vdW interactions were 
taken into account, are shown in Figs[l]-[6j In the Ap- 
pendix, we also show the the PBE and LDA total energy 
curves, Figs [7] - [12] for completion. In all of the Figs 
[T] - [12] the upper panels show the total energy Eq as a 
function of the initial vertical metal atom - graphene dis- 
tance. The total energies were calculated for three ad- 
sorption sites of the metal atom, i.e. top (t), bridge (b), 
and hollow (h) position, respectively. The lower panels in 
Figs [i]- [l2] show the maximal vertical distortion (6 max ) of 
the carbon atoms in the first coordination shell of the ad- 
sorption sites. Positive values of 6 max mean the graphene 
sheet buckles towards the adatom. 

The calculated adsorption energies -E a ds and their cor- 
responding vertical equilibrium distances h eqn are com- 
piled in Tables [I] |III[ and [V] The equilibrium distances 
are the difference between the initial vertical metal atom - 
graphene distance /ij n it and the vertical buckling: h cqu = 
hlnit — b max . Tables [IT] |IV[ and |VI| show the differences 
between the total energies of adsorption sites, e.g. AEj^ b 
between the top and bridge position. These AEo give an 
indication about the preferential diffusion paths of the 
adatoms on graphene. 

Note that all the systems have a total spin moment 
of lfiB due to the single valence electron of the coinage 
metal atoms. This value of the magnetic moment is es- 
sentially preserved even for the shortest distance, show- 
ing that no significant charge transfer takes place between 
the adatom and the graphene sheet. 



We begin by presenting our results for copper on 
graphene. The total energy curves of the PBE+D2 and 
_£.vd\y-DF rjyjj the vdW-DF calculations for this system are shown in 
Figs [I] and [2] (upper panels). For both methods, the top 
position is most preferable. Especially in the PBE+D2 
M]results, the top position competes within a small but no- 
ticeable energy difference with the bridge configuration, 
while the hollow position is definitively not favorable in 
any of the two methods. The PBE (Fig. [7] in the ap- 
pendix) approximation correctly predicts the top site as 
being the ground state, but all the obtained minima are 
comparatively shallow, as one also can see from the ad- 
sorption energies in Table [I] The ordering of the adsorp- 
tion sites obtained with the LDA (Fig. [8]in the appendix) 
are in contrast to the aforementioned vdW and PBE re- 
sults. In LDA, the three sites are in a competitive energy 
range, with the bridge position being the preferred one. 



TABLE I: Adsorption energies (eV) of Cu on graphene, cal- 
culated by four different approximations to the exchange- 
correlation functional, initial vertical Cu-graphene distance h 
and vertical equilibrium distances /i cqu (A), at three binding 
sites, top (t), bridge (b), and hollow (h). 





top 

^ads '^cqu 


bridge 

^ads '^equ 


hollow 

-^ads '''cqu 


PBE+D2 
vdW-DF 
PBE 
LDA 


-0.909 2.11 
-0.684 2.28 
-0.228 2.12 
-0.842 1.96 


-0.901 2.13 
-0.592 2.36 
-0.219 2.08 
-0.874 1.93 


-0.706 2.00 
-0.176 2.57 
-0.063 2.27 
-0.802 1.76 



Table [T] summarizes the adsorption energies -E" a ds and 
corresponding vertical equilibrium distances h cqu for Cu 
and graphene. While PBE+D2 gives E ads of -0.706 to 
—0.909 eV, PBE predicts a much weaker binding, rang- 
ing from —0.063 to —0.228 eV, only. Also, within the 
resolution of 0.1 A PBE and PBE+D2 give similar equi- 
librium distances for the top and the bridge position of 
Cu, but the equilibrium distance at the hollow site is 
12% smaller in PBE+D2 than in PBE. The adsorption 
energies in vdW-DF range from —0.176 to —0.684 eV, 
thereby lying in between those of PBE and PBE+D2. 
Note that this method gives only a very shallow local 
minimum at the hollow site, see Fig. [2] Compared to 
PBE and PBE+D2, the equilibrium distances in vdW- 
DF are shifted towards larger values, exceeding those of 
PBE+D2 by up to 29%. Incidentally, LDA gives binding 
energies quite close to PBE+D2 (within 0.1 eV), but up 
to 12% shorter equilibrium distances. 

In the cases of PBE+D2 and vdW-DF the deforma- 
tion of the graphene sheet is always positive around the 
equilibrium distance, i.e. the carbon atoms and the cop- 
per atom attract each other, see lower panels in Figs [T] 
and [2] For smaller distances, the buckling becomes nega- 
tive, which is obviously a steric effect. Also, the buckling 
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initial vertical Cu - graphene distance (A) 



the actual ground state in LDA. Hence, the energy min- 
imum and buckling do not seem to be correlated in this 
approximation. 



TABLE II: Identifying the diffusion path of Cu on graphene 
from the difference in the total energies A_Bo (meV) at the 
three adsorption sites. 





top to bridge 


top to hollow 


bridge to hollow 




AE t Q - b 




AE b - h 


PBE+D2 


8 


203 


195 


vdW-DF 


92 


508 


416 


PBE 


9 


164 


156 


LDA 


31 


40 


71 



FIG. 1: (Color online) Upper panel: total energy (Eq) as a 
function of the initial vertical Cu adatom - graphene sheet 
distance h for three binding sites, calculated in PBE+D2. 
Lower panel: vertical distortion (& ma x) of the carbon atoms in 
the first coordination shell of the Cu adsorption site. Positive 
values of 6 max mean the carbon atoms buckle towards the 
adsorbed adatom. 
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FIG. 2: (Color online) Upper panel: total energy vs. initial 
vertical Cu - graphene distance for three binding sites, calcu- 
lated in vdW-DF. Lower panel: vertical distortion (6 max ) of 
the carbon atoms closest to the Cu adsorption site. 



is always larger for the top configuration, which is an 
other picture of the increased stabilization of this site 
with respect to the hollow and bridge ones. The verti- 
cal distortions of the graphene sheet obtained from the 
PBE and the PBE+D2 structural relaxations differ by 
less than 10% only. 

The binding order in from the LDA calculations is dif- 
ferent, see Fig. [8j The maximum buckling is obtained 
for a top configuration, although the bridge position is 



From the differences of the total energies at the dif- 
ferent adsorption sites, AEq in Table [TT1 one can obtain 
some insight into the diffusion paths of the coinage metal 
adatoms on graphene, since it has been shown that AEq 
[top - bridge] is the height of the diffusion barrier along 
the carbon-carbon bondsPS Although the energy differ- 
ences vary significantly for AE t Q ^ b in the case of copper, 
the overall trends are identical in PBE+D2 and vdW-DF: 
diffusion will take place along the carbon-carbon bonds, 
since a path over the hollow site is energetically too ex- 
pensive. 



B. Ag on graphene 

The silver adatom on graphene presents quite a differ- 
ent picture than copper and gold do. While PBE+D2 
(Fig. |3| and vdW-DF (Fig. [t| clearly show binding at 
large equilibrium distances, i.e. 2.9 A and 3.3 A, respec- 
tively, pure GGA fails completely to predict any binding, 
see Fig. [9j Note that we repeated the PBE calculations 
with other GGA functional, i.e. PW91 1211 and RPBEP, 
and both failed to predict a minimum as well. As for the 
other two atomic species, LDA gives binding for Ag on 
graphene, with the ordering top-bridge-hollow, see Fig. 
[To} but at least 11% too short equilibrium distances. 

In contrast to Cu and Au on graphene, where one finds 
a small hybridization in the density of states, i.e. a chem- 
ical contribution to the binding, PBE+D2 and vdW-DF 
give a pure physisorption of Ag on grapheneP^ However, 
the energies of the preferred binding sites are ordered 
differently in the two methods. The minima in the total 
energy in PBE+D2 are ordered as hollow - top - bridge, 
while vdW-DF gives top - bridge - hollow. 

The differences in the total energies, listed in Table 
|IV[ of the three binding sites are smaller than 11 meV 
in both vdW methods. In combination with the large 
binding distances, see Table |III[ we conclude that both 
methods predict pure physisorption. Also from Table |TV| 
one can conclude that Ag is the most mobile of the three 
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TABLE III: Adsorption energies (eV) of Ag on graphene, 
calculated by four different approximations to the exchange- 
correlation functional, initial vertical Ag-graphene distance h 
and vertical equilibrium distances h cqu (A), at three binding 
sites, top (t), bridge (b), and hollow (h). 





top 


bridge 


hollow 






rpt it 

ads "'cqu 


-'-'ads '''cqu 


-^ads ''cqu 




PBE+D2 


-0.703 2.90 


-0.700 2.90 


-0.711 2.90 




vdW-DF 


-0.195 3.30 


-0.194 3.40 


-0.192 3.40 


•< 


PBE 


no binding 


no binding 


no binding 


« 

E 


LDA 


-0.368 2.34 


-0.359 2.35 


-0.296 2.59 





•< 




Ag/graphene (PBE+D2) 

* top 
-•- bridge 
— ■- hollow 



2.6 2.8 3.0 3.2 3.4 3.6 
initial vertical Ag - graphene distance (A) 

FIG. 3: (Color online) Upper panel: total energy as a function 
of the initial vertical Ag adatom - graphene sheet distance h 
for three binding sites, calculated in PBE+D2. Lower panel: 
vertical distortion (6 max ) of the carbon atoms in the first co- 
ordination shell of the Ag adsorption site. Positive values 
of femax mean the carbon atoms buckle towards the adsorbed 
adatom. 



coinage metals on graphene and does not have preferred 
diffusion paths as the other two. 

The distortion of the graphene sheet by Ag is negligible 
for the equilibrium distances, see lower panels in Figs 
[3] and [4j supporting the interpretation that it is pure 
physisorption in this system. The buckling away from 
the adatom at smaller distances is again a steric effect, 
as in the other two cases. 

Our PBE+D2 results for Ag on graphene can be com- 
pared to the results in Ref. [25], where van der Waals 
interactions also have been included according to the 
Grimme scheme.^ In contrast to our PBE+D2 results 
in Table III Ref. [23] found smaller adsorption energies, 
i.e. —0.42 to —0.56 eV, and a different binding order. 
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FIG. 4: (Color online) Upper panel: total energy vs. initial 
vertical Ag - graphene distance for three binding sites, calcu- 
lated in vdW-DF. Lower panel: vertical distortion (6 max ) of 
the carbon atoms closest to the Ag adsorption site. 



ences, which might explain the contradicting results be- 
tween our PBE+D2 results and those of Ref. [29] , need 
to be pointed out. First, the value of the Sg coefficient, 
cf. Equ. [4j is not clear from Ref. [33]. Second, the 
Grimme mcthocP^l is parametrized for PBE and not for 
PW91, which has been used in Ref. [29] for the initial 
relaxation. 



TABLE IV: Since Ag is purely physisorpt on graphene the 
difference in the total energies AEo (meV) at the three ad- 
sorption sites are of little significance. 





top to bridge 


top to hollow 


bridge to hollow 




AE t - b 


AEl~ h 


AE b - h 


PBE+D2 


2 


8 


11 


vdW-DF 


1 


3 


2 


PBE 








LDA 


9 


72 


62 



These calculations predicted the top site being more fa- 
vorable than the bridge and hollow site, in accordance 
with our vdW-DF results. Two methodological differ- 



C. Au on graphene 

In the last studied case of a gold adatom on graphene 
all four studied approximations to the exchange- 
correlation energy predict the top adsorption site to be 
most favorable, see upper panels in Figs [5] and [6] as well 
as in Figs [IT] and [12] for comparison. Note that for the 
bridge site vdW-DF predicts a local minimum and PBE 
an inflection point at 2.8 A only. For the hollow site 
vdW-DF gives a shallow minimum at 3.4 A, which is in 
fact lower in energy than the local minimum at the bridge 
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site, while PBE fails to predict binding at this site. Al- 
though, the total energy curve of the PBE+D2 method 
has a minimum at the hollow site as well, the equilibrium 
distance is 3.1 A i.e. even larger than in the case of Ag 
on silver, which indicates physisorption at this site, too. 
Comparing the binding energies of PBE and PBE+D2, 
the latter essentially adds —0.8 eV to E a( ^ B of the former. 
This indicates the importance of the vdW interaction to 
the binding, which is even more pronounced in this sys- 
tem than for Cu on graphene. 



TABLE V: Adsorption energies (eV) of Au on graphene, cal- 
culated by four different approximations to the exchange- 
correlation functional, initial vertical Au-graphene distance h 
and vertical equilibrium distances /i oq u (A), at three binding 
sites, top (t), bridge (b), and hollow (h). Note that vdW-DF 
and PBE do not predict binding of Au at the hollow site. 





top 

-^ads '^cqu 


bridge 


hollow 

rph j h 
-^ads '*equ 


PBE+D2 
vdW-DF 
PBE 
LDA 


-0.886 2.53 
-0.385 2.65 
-0.099 2.54 
-0.732 2.22 


-0.881 2.79 
-0.314 2.72 
-0.081 2.72 
-0.698 2.26 


-0.870 3.10 
-0.322 3.40 
no binding 
-0.451 2.40 
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FIG. 5: (Color online) Upper panel: total energy as a function 
of the initial vertical Au adatom - graphene sheet distance h 
for three binding sites, calculated in PBE+D2. Lower panel: 
vertical distortion (6 max ) of the carbon atoms in the first co- 
ordination shell of the Au adsorption site. Positive values 
of fcmax mean the carbon atoms buckle towards the adsorbed 
adatom. 
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For the calculated adsorption energies, Table [V] we ob- 
tained trends that are similar to those of Cu on graphene, 
cf. Table [I] The PBE+D2 method predicts E ads ranging 
from -0.870 to -0.886 eV, which exceed even the LDA 
results for this system. PBE, on the other hand, pre- 
dicts significantly lower -E a d s of less than -0.099 eV, only, 
while vdW-DF gives adsorption energies that lie in be- 
tween those of PBE and PBE+D2, ranging from —0.310 
to -0.385 eV. 



The equilibrium distances h cqu shown in Table 



V] agree 

rather well for PBE+D2 and vdW-DF, i.e. withdevia- 
tions of less than 9%. Note that the buckling of the 
graphene sheet is 50% smaller in PBE+D2 than in vdW- 
DF, i.e. PBE. Even PBE gives comparable h eqn for the 
top and bridge site, while LDA overbinds and predicts 
h cqu that are up to 23% shorter than those obtained in 
PBE+D2. 



TABLE VI: Identifying the diffusion path of Au on graphene 
from the difference in the total energies AEo (meV) at the 
three adsorption sites. 



Au/graphene (vdW-DF) 

* top 
-•- bridge 
— ■- hollow 




2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 
initial vertical Au - graphene distance (A) 

FIG. 6: (Color online) Upper panel: total energy vs. initial 
vertical Au - graphene distance for three binding sites, calcu- 
lated in vdW-DF. Lower panel: vertical distortion (6 max ) of 
the carbon atoms closest to the Au adsorption site. 





top to bridge 


top to hollow 


bridge to hollow 




AEl- b 


AE f - h 


AE b - h 


PBE+D2 


5 


16 


11 


vdW-DF 


71 


63 


8 


PBE 


18 






LDA 


34 


281 


247 



As in the case of copper, the differences of the total 
energies between the top and bridge site vary greatly 
among the employed vdW approximations. Still, one can 
conclude that the gold adatom is likely to diffuse along 
the carbon-carbon bonds. Therefore, the usage of PBE 
is justified for studies of, for instance, the mobility and 
clustering of gold on graphene as in Ref . [35J . 
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IV. SUMMARY AND CONCLUSIONS 

We performed a systematic DFT investigation of 
the adsorption of copper, silver, and gold adatoms on 
graphene, especially taking van der Waals interactions 
by the vdW-DF and the DFT+D2 methods into ac- 
count. For copper and gold we found that the PBE 
parametrization to the exchange-correlation energy pre- 
dicts the same ordering of the adsorption sites as vdW- 
DF and PBE+D2, i.e. it also favors the top over the 
bridge and hollow positions. We also find that the 
non-local interactions increase the calculated adsorption 
energies of Cu and Au on graphene by up to 0.8 eV 
(PBE+D2) and 0.45 eV (vdW-DF), respectively. 

The predicted vertical equilibrium distance calculated 
with PBE, for Cu and Au adsorption, agrees to bet- 
ter than 13% with the more advanced non-local meth- 
ods. Taking vdW interactions into account during the 
structural relaxations in the PBE+D2 method did not 
significantly change the buckling of the graphene sheet 
in the case of Cu and Ag compared to PBE calcula- 
tions. Compared to the PBE calculation, we found a 
50% smaller buckling of the graphene sheet in the case 
of gold with the PBE+D2 method that might be related 
to its parametrization. Ripples in the graphene structure 
give rise to both effective scalar and vector potentials in 
the electronic structure of graphene, see e.g Ref. [53] . 
The adsorption of Cu or Au atoms may be an important 
tool to enhance this rippling and hence the strength of 
these scalar and vector potentials. For the adsorption 
of silver on graphene our calculations suggest that it is 
purely of van der Waals type. Different generalized gra- 
dient approximations to the exchange-correlation func- 
tional fail to give any binding at all, while the PBE+D2 
and vdW-DF predict physisorption at large equilibrium 
distances of 2.9 and 3.3 A, respectively. Also the calcu- 
lated distortion of the graphene sheet upon adsorption of 
a silver atom was found to be negligible. From the differ- 
ences in the total energies we conclude that diffusion of 
Cu and Au takes place along the carbon-carbon bonds, 
while the Ag adatoms can diffuse almost unrestricted on 
the graphene sheet. 
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V. APPENDIX 

In the appendix we present calculated total energy 
curves and vertical distortions of the carbon atoms clos- 
est to the adsorption sites for the PBE and LDA func- 
tional. The Cu results are found in section A, the Ag 
results in section B and the Au results in section C. 



A. Cu/graphene (PBE and LDA) 




1.6 2.0 2.4 2.8 3.2 3.6 
initial vertical Cu - graphene distance (A) 



FIG. 7: (Color online) Upper panel: total energy vs. initial 
vertical Cu - graphene distance for three binding sites, calcu- 
lated in PBE. Lower panel: vertical distortion of the carbon 
atoms closest to the Cu adsorption site. 
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1 .6 2.0 2.4 2.8 3.2 3.6 
initial vertical Cu - graphene distance (A) 



FIG. 8: (Color online) Upper panel: total energy vs. initial 
vertical Cu - graphene distance for three binding sites, calcu- 
lated in LDA. Lower panel: vertical distortion of the carbon 
atoms closest to the Cu adsorption site. 



B. Ag/graphene (PBE and LDA) 



C. Au/graphene (PBE and LDA) 



-460.8 

_ -461.0 
> 

-461.4 
0.0 

^ -0.2 

CD 

E 
£} 

-0.4 



■ 1 1 1 


1 1 1 1 1 
^°°*~*" " ■ ■ ■ ■ ■ 


-i | 1 1 1 1 1 1 , p- 


1 , 1 


Ag/graphene (PBE) 

* top 
— •- bridge 

■ hollow 
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FIG. 9: (Color online) (Color online) Upper panel: total en- 
ergy vs. initial vertical Ag - graphene distance for three bind- 
ing sites, calculated in PBE. Lower panel: vertical distortion 
of the carbon atoms closest to the Ag adsorption site. 
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2.0 2.4 2.8 3.2 3.6 

initial vertical Au - graphene distance (A) 

FIG. 11: (Color online) Upper panel: total energy vs. initial 
vertical Au - graphene distance for three binding sites, calcu- 
lated in PBE. Lower panel: vertical distortion of the carbon 
atoms closest to the Au adsorption site. 
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Ag/graphene (LDA) 




* top 
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1,1,1,1 



2.0 2.4 2.8 3.2 3.6 

initial vertical Ag - graphene distance (A) 

FIG. 10: (Color online) (Color online) Upper panel: total 
energy vs. initial vertical Ag - graphene distance for three 
binding sites, calculated in LDA. Lower panel: vertical dis- 
tortion of the carbon atoms closest to the Ag adsorption site. 
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FIG. 12: (Color online) Upper panel: total energy vs. initial 
vertical Au - graphene distance for three binding sites, calcu- 
lated in LDA. Lower panel: vertical distortion of the carbon 
atoms closest to the Au adsorption site. 
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